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ABSTRACT 

A sample consisting of 27 X-ray selected galaxy clusters from the XMM-LSS survey is used 
to study the evolution in the X-ray surface brightness profiles of the hot intracluster plasma. 
These systems are mostly groups and poor clusters, with temperatures 0.6-4.8 keV, spanning 
the redshift range 0.05 to 1.05. Comparing the profiles with a standard /3-model motivated 
by studies of low redshift groups, we find 54% of our systems to possess a central excess, 
which we identify with a cuspy cool core. Fitting /3-modcl profiles, allowing for blurring by 
the XMM point spread function, we investigate trends with both temperature and redshift in 
the outer slope (j3) of the X-ray surface brightness, and in the incidence of cuspy cores. Fits to 
individual cluster profiles and to profiles stacked in bands of redshift and temperature indicate 
that the incidence of cuspy cores does not decline at high redshifts, as has been reported in 
rich clusters. Rather such cores become more prominent with increasing redshift. j3 shows 
a positive correlation with both redshift and temperature. Given the (3-T trend seen in local 
systems, we assume that temperature is the primary driver for this trend. Our results then 
demonstrate that this correlation is still present at z > 0.3, where most of our clusters reside. 

Key words: galaxies: clusters: general - galaxies: X-ray. 



1 INTRODUCTION 

Clusters of galaxies, as the largest virialised gravitationally-bound 
products of the process of hierarchical structure formation, are 
powerful probes for both testing cosmological models and tracing 
structural evolution (e.g., Voit 2005). One of the most important 
properties of galaxy clusters is their mass. Since cluster mass can- 
not be directly observed, it is studied indirectly through observables 
such as X-ray radiation emitted by the intracluster medium (ICM) 
which represents 80% of the total baryonic component of galaxy 
clusters at z=0 (Ettori et al. 2004) and accounts fo r about 10% o f 
the total (including dark) mass content of clusters (Sarazin 1986). 
The study of the ICM can provide important insights into the evo- 
lution and dynamics of cluster and their member galaxies. 

Observationally, there are two distinct classes of clusters: cool 
core (CC) clusters with dense gaseous core regions in which gas 
temperature drops inwards, and non-cool core (NCC) clusters with 
shallower core profiles which o ften exhibit more internal struc- 
ture ( e.g., [jpnes & Forman 1984; Ota & Mitsudal l2004; Per es et all 
1 19981 and JSchuecker et alj|200lh . Cool core clusters have sharply 
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peaked X-ray emission at their centres due to the rise in central gas 
density which accompanies the central cooling. However, the gas 
is not observed to cool to very low temperatures at the rates naively 
expected from the observed core X-ray luminosities, and a con- 
sensus has now emerged that this is due to the effects of feedback 
from a central active galactic nucleus, which limits the effects of 
cooling through processes which are still not very we ll understood . 
For reviews of cool cores in clusters, see for example. iFabianl 19941 
iDonahue & Voidl2004l and lPeterson & Fabianll2006l . 

In the local Universe, some stud ies have found that nearly two 
thirds of clusters have cool core s (e.g. lPeres et alJ 1998. Whi te et al.l 
1 1 9971 andl Vikhlinin et al. 2007). Ho wever, other studies gave differ- 
ent values: Edge. Stewart & Fabian] 1 19921) found a CC fraction as 
high as 90%, while the results of IChen et all d2007t> indicated that 
49% of local clusters host cool cores. These differences relate to 
both the selection of the cluster sample, and the way in which cool 
cores are identified within them. 

The evolution of cooling within clus ters provid es an important 
probe of the history of cosmic feedback (Voit 2005). At intermedi- 



ate redshifts (z ~ 0.15 — 0.4). lBauer et al.l 12005). using a sample 
of 38 X-ray-luminous clusters, found that cool cores appeared still 
to be common, with an incidence nearly identical to that in lumi- 
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nous low-redshift clusters. Consequently, they suggested that heat- 
ing and cooling processes must have stabilised in massive clusters 
since z ~ 0.4. 

At higher redshifts, Ivikhlinin et al.1 J2007I) reported that the 
fraction of clusters with cuspy X-ray cores drop ped from ~70% 
at 2 ~ to ~15% at z > 0.5. ISantos et alj J200§) compared the 
fraction of clusters with non-cool cores, moderate cool cores and 
strong cool cores in nearby (0.143 < z < 0.3) and high redshift 
(0.7 < z < 1.4) clusters. These authors detected a significant 
fraction of clusters harbouring moderate cool cores out to z=1.4, 
similar to the fraction in their low-redshift sample. However, they 
noticed an absence of clusters with strong-cool cores at redshift 
z > 0.7. 

Regarding the spa tial distribution of the ICM, 

ICavaliere & Fusco-Femiand < 19761) introduced the /3-model 
profile, motivated by the distribution expected for an isothermal 
plasma in hydrostatic equilibrium with a virialised mass distribu- 
tion. Although it is now known that the gas is rarely isothermal, the 
/3-model is generally found to give a reasonable representation o f 
X-ray surface brightness profiles (e.g. iNeumann & Ar naud 1999). 
However, an additional central component is usually required to fit 
the inner regions of CC clusters (e.g. IPratt & Ar naud 2002), and 
detailed studies of surface brightness profiles of clusters extending 
to large radii have shown that the logarithmic slope con tinues to in- 
creas e slowly towards large r rad ii (e.g . Vikhlinin, Forman & Jonesl 
1 19991 ICroston et all 120081 and IMaughan et alj |2008h . It is less 
clear whether this prog ressive steepening is also prese nt in galaxy 
groups. For example, iRasmussen & Ponmanl J2004D traced the 
surface brightness out to R500 in two rich groups and found them 
to be well fitted by simple /3-models the whole way. 

For /3-model fits, the j3 parameter, which characterises the 
outer slope of th e surface brightness pro file, has a value of « 2/3 
for rich clusters Uones & Formanll 19841) and lower values for poor 
cluste r s and galaxy groups {Fi noguenov, Reiprich & Bo hringen 
I2001L iHelsd on & Ponman 2000 and Ho rner et all |l999l) . Sev- 
eral studies have shown that f3 has a mild positive trend 
with the average temperature of th e ICM in nearby clus- 
ters i Vikhlinin Forman & Jonesl 1 19991 ICroston et all 120081 and 
IMaughan et alj|2008l) and that the value in poorer galaxy groups 
is lower (i.e. flatter surfac e brightness slope) than that in clusters 
JOsmond & Ponmanl2004l) . In terms of evolution, a study of Chan- 
dra data for 115 clus ters spanning the range 0.1 < 2 < 1.3 by 
IMaughan et alj ([2008) shows some indication that the /3-T correla- 
tion is weaker for clusters at 2 > 0.5. 

Low-mass galaxy clusters or groups, with ICM tempera- 
tures less than 2-3 keV, play an important role in the evolu- 
tion of galactic systems because they lie at a transition be- 
tween the field environment and rich cluster environments, and 
also because non-gravitational processes have a larg er impact in 
groups than in rich clusters (e.g. , Zabludoff & Mulchae v 1 19981 
IPonman. Sanderson & Finoguenovl 120031 and ISun et alj 2009). 
However, these poor systems, and in particular the evolution of 
their properties, have received rather little attention. This is mainly 
due to the difficulty in detecting and studying them, especially at 
large redshifts, due to their faint X-ray emission and small comple- 
ment of galaxies. 

For these reasons, research on the evolution of galactic sys- 
tems in the regime of groups and poor clusters has only started re- 
cently, as a result of improvements in observing capabilities in both 
the X-ray and optical. By comparing optically-selected systems at 
0.3 < z < 0.55 with nearby groups. Iwilman et al.1 J20 05) showed 
that the fraction of group members undergoing significant star for- 



mation increases strongly with redshi ft out to z ~ 0.45. H owever, 
the study of X-ray selected groups bv ljeltema et alj J2007I) showed 
a contrary result: they did not observe significant evolution in the 
morphology or star formation of the galaxy populations in their 
0.2 < z < 0.6 groups compared to low-redshift X-ray luminous 
groups. They argued that this discrepancy could be due to differ- 
ent selection methods, since optically-selected systems are typi- 
cally lower in mass and contain more spiral galaxies and therefore 
a stronger evolution in the galaxies is expected. They also found 
that their moderate redshift groups had galaxy populations simi- 
lar to clusters at the same redshift; in particular, a large fraction of 
early-type galaxies and a low fraction of galaxies with significant 
star formation. However, in contrast to the situation in low red- 
shift X-ray bright groups, a sig nificant fraction of these intermedi- 
ate re dshift groups were found iMulchaev et al. 2006, J eltema et al.l 
1 20071) to have no bright early-type galaxy at the centre of the X-ray 
emission, or to have a central galaxy with multiple nuclei. This was 
taken as evidence for the dynamical youth of many of these groups. 

The small number of studies which have addressed the evo- 
lution in the X-ray properti es of galaxy groups have found little 
convincing evidence for any. Jelt ema et alj d2006l) . in a multiwave- 
length study of six galaxy groups and poor clusters at intermedi- 
ate redshift (0.2-0.6), found that they appear to follow the scal- 
ing relations between luminosity, temperature, and velocity dis- 
persio n defined by low-r edshift groups and clusters. This is also 
true ( Jeltema et alj 120091) for three higher redshift poor clusters 
from the AEGIS survey. A study of evolution in the L-T relation 
based on the present XMM-LSS cluster sample by IPacaud et al.1 
(2007), taking into account the selection function of the survey, 
found that the range of models consistent with the data included 
self-similar evolution, an d also (marginally) a no-evolution model. 
iFinoguenov et alj l2007l) extracted a larger sample of 72 clusters 
(mostly poor ones) from the XMM-Newton observations of the 
COSMOS field, and found no evidence for evolution in the lumi- 
nosity function of these systems out to 2 ~ 1, though the quality 
of their data did not permit them to study the morphology of the 
X-ray emission. 

Motivated by the the paucity of information available for the 
evolution of the ICM in the important environment of low-mass 
galaxy clusters, we aim in this study to shed light on the spatial 
distribution of the ICM in X-ray selected clusters covering a wide 
redshift range (2 ~ — 1), paying special attention to trends in the 
slope and central cuspiness of the X-ray emission. 

The paper is constructed as follows: in section 3.2, we describe 
the data and briefly introduce the properties of the cluster sample; 
then we describe the data reduction used to produce X-ray surface 
brightness profiles. In section 3.3, we present our results, starting 
with the individual cluster profiles, and then profiles of redshift- 
stacked and temperature-stacked clusters. In section 3.4, we discuss 
the implication of our results and compare them with other studies. 
Finally, in section 3.5, we summarise our conclusions. 

Throughout this paper, we adopt the co smological parameter s 
from the five-year WMAP data presented bv lHinshaw et al.1 (2009), 
namely: H = 70.5 km s _1 Mpc -1 , Q. m = 0.27, fi A = 0.73. 



2 DATA 

2.1 The sample 

Our sample is based on the 29 Class 1 (CI) clusters from the X- 
ray Multi-Mirror Large-Scale Structure (XMM-LSS) survey. The 
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XMM-LSS CI cluster sample is a well-controlled X-ray selected 
and spectroscopically confirmed cluster sample. The criteria used 
to select the members of this sample guarantee negligible contam- 
ination of point-like sources. The observations of the clusters were 
performed in a homogeneous way (10-20 ks exp osures). For full 
details of the CI sample, see lPacaud et alj ( 120071) . Detailed infor- 
mation on the selec tion function of the CI sample can be found in 
IPacaud et all feOOd) . 

The C 1 sample is dominated by groups and poor clusters with 
temperatures of 0.63 < kT < 4.80 keV, spanning a redshift range 
0.05 < z < 1.05. Typically, we have a few hundreds X-ray counts 
for each cluster, with only a few having over a thousand detected 
photons. Two of the 29 clusters with less than 80 counts had to be 
excluded from our analysis because their data were inadequate for 
useful profiles to be extracted. The excluded clusters are XLSSC 
clusters 39 and 48. Hence our sample consists of 27 clusters. Clus- 
ter 47, with 81 counts, was a marginal case. We were unable to 
constrain a fit to its individual profile, but its data were included in 
the analysis of the stacked profiles. Key properties of the sample 
are presented in Table[T] 

2.2 Data analysis 

To construct the X-ray surface brightness profiles for each clus- 
ter, we used three X-ray FITS images, three exposure maps and 
one segment ation map, all produc ed using the production pipe line 
described in IPacaud et alj ([2006). The images were taken by the 
MOS1, MOS2 and PN imagers on board the XMM-Newton satel- 
lite in the energy band [0.5-2.0] keV with exposure times ranging 
from 10 to 20 ks. The exposure maps are FITS images containing 
the vignetting-corrected exposure of the clusters as a function of the 
sky position. A single segmentation map, generated by SEXTRAC- 
TOR was used for each cluster to remove contaminating sources. 

The right ascension (RA) and declination (Dec) values of the 
centres were determined as outlined in Pacaud et al. (2006). But 
when we examined the X-ray profiles of the clusters, some showed 
dips at the centre. For these clusters, we mosaicked the three im- 
ages, smoothed the resulting image and took the coordinates of the 
pixel with the maximum photon counts and modified the cluster 
centre accordingly. The modified centres at most are only ~14 arc- 
seconds from the original values but remove the central dips in the 
profiles. The clusters with modified centres have XLSSC numbers: 
50,28,40,1,47 and 5. The RAs and Decs in Tableware the modi- 
fied centres and the original coordinates can be found in Table 1 in 
IPacaud et all fc007l) . 

Since the angular size of our clusters is small, background 
removal using a local estimate works well. The background was 
taken from an annulus extending from 2xRsoo to 3xRsoo about 
each cluster, where R500 is the radius within which the mean clus- 
ter mass density is 500 tim es the critical densit y of the Universe 
at the cluster redshift. As in lPacaud et alj 120071) . the R500 values 
were calculated using 

Rsoo = 0.388 x T ' 63 x h 70 . 5 (z)- 1 Mpc, (1) 

where T is the ICM temperature in keV and hro.5 is the 
Hubble constant in units of 70.5 kms - Mpc - . This 
formula was originally deriv ed from M-T relation of 
iFinoguenov. Reiprich & Bohring er (2001). 

Two background components were evaluated: the photon 
background component and the particle background component. 
These were separated using the fact that photons are vignetted, 
whilst particles are not. Hence the relationship between count rate 



and effective area for pixels in the background annulus gives an es- 
timate of the photon background component from the slope, and of 
the particle background component from the intercept. 

Surface brightness profiles were extracted from a standard set 
of annuli (as a fraction of R500 to facilitate later stacking), extend- 
ing to 3xRsoo- For each annulus we removed the particle back- 
ground and computed the the vignetting-corrected count rate (in 
ct/s/pix) and its error. The MOS1, MOS2 and PN profiles, gen- 
erated in this way were then combined, and their errors added in 
quadrature. Profiles were extracted up to 3 x R500 where they flat- 
tened and reached the photon background values. The final column 
in Table[TJis the total X-ray counts within 3 x R500 after subtracting 
from it both the photon and the particle background. 

The X-ray cameras on the XMM-Newton satellit e have a point 
spread fun ction (PSF) of ~ 6 arcsec FWHM, see IStriider et alj 
J2001I) and lTurner eTatl J200lh . Correction for PSF blurring is im- 
portant to avoid biased estimation of the parameters of the cluster's 
radial profile, especially, th e core radius, r core . W e applied the PSF 
correction method used by lArnaud et alj J2002f) (and described in 
detail therein) which analytically computes a photon redistribution 
matrix (RDM) based on the properties of the three XMM-Newton 
cameras and depends on the energy band used and off-axis angle 
between the centre of the camera and the cluster position. The PSF 
matrices for the three cameras were weighted by the source counts 
in each camera and combined to produce a matrix appropriate to 
the coadded profiles. 

We fitted the surface bright ness profile with a /3-model 
dCavaliere & Fusco-Femianoll 19761) 



S(r) = 5o(l + (r/r corc ) 2 ) 



-3/3+0.5 



(2) 



where So is central brightness (cts/s/pix) and r corc is the core radius 
(in units of Rsoo)- The model was blurred with the PSF redistribu- 
tion matrix, and fitted to the surface brightness data. The best val- 
ues of r corc and were estimated by computing the minimum x 2 
value on an adaptively refined r cola -/3 grid, la errors were com- 
puted for r coro and /3, and la, 2a and 3a error regions computed 
in the r corc -/3 plane for each fit. 

As will be seen below, a number of systems show a central 
excess above the fitted /3-model. Such a central cusp suggests the 
possible presence of a cool core. Since a central excess may distort 
the /3-model fit, we attempted to fit a model with the central bin 
excluded, but given the limited statistical quality of our data, loss 
of the central bin resulted in poorly behaved fits in many cases. 
We therefore adopted the approach of fitting a /3-model with core 
radius fixed at a value (as a fraction of R500) motivated by the ob- 
served profiles of local groups in which detailed modelling of the 
surface brightness has been possible. A central excess above this 
model then indicates the presence of a cuspy core. 

The fixed value of the core radi us we adopt as a canonica l 
value for poor clusters is taken from iHelsdon & Ponmanl (2000), 
who studied 24 X-ray-bright galaxy groups. For half of their sam- 
ple they found that two-component /3-models were required to give 
acceptable fits to the surface brightness distribution. The outer com- 
ponent represented the intragroup gas, whilst core emission could 
be distinguished by a clear shoulder in the profile in many cases, 
and was fitted by the inner /3-model component. The median value 
of r CO re for the outer component in the 12 clusters was found to 
be 60 kpc. Correcting to our va lue of Hp, this median r cor e would 
be 42.6 kpc. The 12 systems in lHelsdon & Ponmanl ( 120000 had an 
average temperature of 1.07 keV. To calculate this r COTC as a frac- 
tion of -R500, we used the R,5pp(T) equation above, which gives 
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XLSSC R.A. Dec. Redshift kT r 5 oo r corc /i? 5 oo Central Excess Counts 

number (J2000) (J2000) (keV) (Mpc) Fitted r core Fixed r corc Factor (/ c ) 
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Table 1. List of the properties of the 29 C 1 galaxy clusters sample sorted according to their redshifts ffacaud et al]|2007l) with the fitted parameters of the 26 
successfully constrained clusters. The and r cor e values in the seventh and eighth columns are for the the /3-model with both r cor e and freely fitted. The 
values in the ninth column are for the /3-model with fixed r core =0.105 X Rsoo- All errors are 1<t errors. Clusters 47 failed to converge but its data were 
not excluded from the stacked analysis. The central excess factor, f c (tenth column), is calculated by dividing the observed surface brightness by the predicted 
surface brightness (from the fitted model) for the innermost radial bin, that is 0.05 X R500 • If fc is greater than 1 then this is considered as an indication that 
the cluster has a cool core and vice versa. The stars (***) denote unconstrained errors. 



r core =0.105 x Rsoo- This value of r core was therefore used in our 
fixed-core fits. 

To quantitatively determine whether a profile of a cluster (or 
a stacked set of clusters) has a central brightness excess, and there- 
fore a CC, we define the central excess factor (/ c ) as the ratio be- 
tween the observed surface brightness and predicted surface bright- 
ness (from the fitted fixed-core model) for the innermost radial bin, 
at r = 0.05 x Rsoo- If fc is greater than 1 then this was considered 
as an indication that the cluster has a CC and vice versa. The er- 
ror on f c is derived simply from the error on the innermost surface 
brightness value. We checked the effects of excluding the central 
bin for the fixed-core fits. This has little effect on the fitted /3 value, 
but in cases where there is a central excess, it results in a somewhat 
lower normalisation for the fitted model, and hence a slightly (up 
to 10%) higher value of f c . 



2.3 Stacked profiles 

The statistical quality of our individual profiles is limited, so we 
stacked the observed profiles of clusters with similar redshifts and 
temperatures, producing higher quality profiles which might high- 
light any trends with temperature or redshift. 

Since each cluster has a different R500, we extracted the pro- 



files in fixed radial bins in units of R500 up to 3 x Rsoo- The range 
from (O-l)xRsoo was divided into 20 equally spaced bins, from 
(l-2)xR5oo into 15 equally spaced bins and the range from (2- 
3)xRsoo into 10 equally spaced bins. These different bin widths 
allow for the decline in flux with radius whilst keeping the inner 
bins sufficiently fine to resolve the core. This distribution of bins 
was chosen after some experimentation to obtain the best fit con- 
straints. 

Before the stacking, the profiles were multiplied by the scale 
factor: 

I (3 ) 

where 

A = R500 (4) 

to account for the cluster line of sight depth, 

to correct for the change in critical density of the Universe and 
C={l + zY A (6) 

to eliminate the effect of cosmological dimming. The aim of the 
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scaling is to allow for the effects of variable cluster depth and for 
cosmological factors, so that all profiles would be similar (to within 
the rather weak temperature dependence of the X-ray emissivity) 
in the case where clusters are simple self-similar systems, evolving 
with the critical density of the Universe. 

The profiles for each component cluster were then added bin 
by bin and their errors quadratically summed to generate a stacked 
profile. The photon background values for each cluster were scaled 
by the same factors as the source profiles, before being combined. 
This coadded photon background was then included in the fitted 
model as a fixed background level. The PSF redistribution matrices 
were weighted by the scaled count rate for each cluster before being 
combined to produce a composite matrix. 



the intermediate-^ stacked set, with average redshift and tempera- 
ture of 0.30 and 1.69 keV respectively. The high-z set contains 9 
clusters spanning a redshift range 0.43 to 1.05 (average 0.72) and 
having temperatures from 2.20 to 4.80 keV (average 3.51 keV). 
The profiles and the error contours of the three stacked sets are 
presented in Fig. [4] and the fitted parameters for the free and fixed 
/3-model fits with la errors are shown in the first three rows of 
Tabled 

The central excess factor, f c is seen to increase with redshift; 
for the low-z stack it is 1.30, increasing to 1.56 and 1.95 for the 
intermediate and high z systems respectively. Table [2] also shows 
that /3 (for both free and fixed r CO rc fits) increases with redshift, 
whilst the r corc values are essentially constant. 



3 RESULTS 

3.1 X-ray surface brightness profiles of individual CI 
clusters 

Most individual CI clusters profiles fit successfully with a free- 
r core /3-model, and nearly all have well-constrained fixed-r corc /3- 
model fits. The profile of cluster 47, with only 81 counts, is the only 
one for which we could not achieve a useful /3-model fit. For clus- 
ters 51 (160 counts), 49 (157 counts) and 5 (130 counts), although 
best fit models were obtained, the upper bound of the free-r- corc j3 
values were not constrained, and no r COIC -/3 contour plots could be 
produced. 

The surface brightness profiles with fits and associated la, 2a 
and 3cr error contours for individual clusters with redshift ranges 
0.05-0.17, 0.26-0.33 and 0.43-1.05 are shown in Fig.[j][2]and[3]re- 
spectively, and numerical results are given in Table Q] The profile 
of cluster 50 is unusual; it has a remarkably large r corc value, of 
1.12xRsoo (see Fig.[TJ, and an elongated spatial extension indicat- 
ing a cluster in a state of merging, with a highly unrelaxed core. 

Amongst the 26 clusters with constrained fixed-r corc fits, 21 
(81%) possess CCs according to the criterion outlined above (i.e. 
f c > 1), whilst the remaining five (19%) are non-cool core (NCC) 
(f c < 1) systems, as shown in TableQ] However, for some of these, 
the classification must be regarded as uncertain, since the error bar 
on fc crosses unity. This is the case for 7 (from the 21) CC systems, 
and 2 (of the five) NCC clusters. Hence, at the la level, 54% (14 
out of 26) of our clusters show a central excess in surface brightness 
which may indicate the presence of a cool core. 

The median value of /? for free-r CO rc fits to the 26 clusters 
is 0.61. A similar median value (/3=0.63) is obtained from the 
fixed-r C orc fits. As for the r corc , its median value is 0.08xR,5oo- 
As expected, this is rather smaller than the canonical value of 
0.105 xRsoo for the group scale component of the emission, due 
to the influence on the fits of cuspy cores in many systems. 



3.2 X-ray surface brightness profiles of redshift-stacked 
clusters 

The C 1 clusters span a redshift range of 0.05 to 1 .05. To probe how 
the X-ray surface brightness profiles evolve with redshift, we di- 
vided the CI sample into three redshift ranges: 0.05-0.17 (low-z), 
0.26-0.33 (intermediate-z) and 0.43-1.05 (high-z). The low-z set 
consists of 6 clusters with an average redshift of 0.11 and a tem- 
perature range from 0.63 to 3.50 keV (average 1.33 keV). Only 
one of these (cluster 50) has kT > 2 keV. Twelve clusters fall in 



3.3 X-ray surface brightness profiles of temperature-stacked 
clusters 

The full temperature range of our C 1 sample (0.63 to 4.80 keV) was 
divided into three subsets. The coolest set (0.63 keV < kT < 1.34 
keV) contains ten clusters with average fcT=1.06 keV and average 
2=0.20. There are nine clusters in the second set ( 1 .60 ke V < kT < 
2.80 keV) with averages kT=2.13 keV and 2=0.38, and the hottest 
set (3.20 keV < kT < 4.80 keV) contains eight clusters with av- 
erage temperature and redshift of 3.76 keV and 0.68 respectively. 
The stacked profiles and the associated la, 2a and 3a contours are 
shown in Fig.|5]and the fitted parameter values in Table|2] 

All three temperature-stacked sets show evidence for CCs, 
with f c > 1. However, f c does not show a monotonic trend with 
temperature as was the case for the redshift-stacked clusters. The 
intermediate-temperature set shows the strongest central excess, 
with / c =2.17±0.28. Similarly, the /3 values, for both fixed and 
free r core fits, do not show a monotonic trend across the full T 
range of our sample, although it is clear that the hotter systems 
(kT > 1.5 keV) have /3 values significantly higher than the groups 
in our coolest bin. The core radius, r corc appears remarkably stable 
in these stacked clusters, fitting at a value 0.07xRsoo. This is also 
essentially the case for the redshift-stacked clusters, in which the 
high and low redshift stackes fitted at r corc = 0.07xRsoo whilst 
the intermediate-redshift stack gives r corc =0.06xR5oo, which is 
the same within errors. 



4 DISCUSSION 

4.1 X-ray surface brightness profiles of z ~ 0.3 clusters 

IPacaud et al.l (2007J) show (see their Fig. 3) that the redshift distri- 
bution of the CI clusters, which spans the redshift range 0.05 to 
1.05, has a pronounced peak around z ~ 0.3. More than 40% (12 
systems) of our clusters are concentrated in the relatively small red- 
shift range 0.26 < z < 0.33. The average temperature of these 12 
clusters is 1.69 keV. Their average M 50 o is 3.96 xlO 13 M . This 
puts them in the realm of groups or poor clusters. To our knowl- 
edge, this is the best sample of X-ray selected groups at 2 ~ 0.3 
studied to date, and hence our individual and stacked X-ray pro- 
files of these clusters provide the best available X-ray profile of 
low-mass clusters at intermediate redshift, and should be useful for 
future comparative studies. The individual X-ray profiles with the 
the r corc -/3 contours for these cluster are shown in Fig. [2] and the 
stacked profile is the second panel in Fig. [4] 

The free-r cor c fit to the stacked profile for these clus- 
ters gives p = 0.51 (see Table [2] second row) which is 
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Figure 1. X-ray surface brightness profiles of the individual CI clusters with redshift 0.05 < z < 0.17, ordered according to redshift and the associated 
ltr, 2f7 and 3<r contours. The dashed lines are the fitted /3-model profiles with both r CO rc and (3 freely fitted, while the solid lines are for the fitted profiles 
with free (3 and r C ore fixed to 0.105 x R500- 
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Figure 2. X-ray surface brightness profiles of the individual CI clusters with redshift 0.26 < z < 0.33, ordered according to redshift and the associated 
constrained la, 2cr and 3ct contours. The dashed lines are the fitted /3-model profiles with both r cor c and freely fitted, while the solid lines are for the fitted 
profiles with free /3 and r cor c fixed to 0.105 X R500. 



8 A. Alshino et al. 




Figure 3. X-ray surface brightness profiles of the individual CI clusters with redshift 0.43 < z < 1.05, ordered according to redshift and the associated 
constrained lcr, 2a and 3<r contours. The dashed lines are the fitted /3-model profiles with both r CO rc and freely fitted, while the solid lines are for the fitted 
profiles with free and r cor c fixed to 0.105 X Rsoo- 
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Range 



Stacked Clusters 
XLSSC 



Average z 



Average kT 
(keV) 



Fitted r cor e 



f'corc/-R500 





Fixed T-corc 



Central Excess 
Factor (/ c ) 



Redshift-stacked clusters (all clusters) 



z:0.05-0.17 11,52,21,41,50,35 0.11 1.33 0.46 

z:0.26-0.33 25,44,51,22,27,8,28,13,18,40,10,23 0.30 1.69 0.51 

z:0.43-1.05 6,36,49,1,2,47,3,5,29 0.72 3.51 0.60 



+0.01 
0.01 

+ 0.01 
0.01 

+0.06 
0.04 



0.07 
0.06 
0.07 



+0.01 
0.01 

+0.01 
0.01 

+0.02 
0.02 



W.003 
-0.006 



0.49Z 
5Q0+ 0009 

U.O9U_ Q11 

0.70"* 



hO.02 
-0.03 



1.30 ±0.09 
1.56 ±0.11 
1.95 ±0.25 



Temperature-stacked clusters (all clusters) 



kT: 
kT: 
kT: 



0.63-1.34 52,11,21,13,35,51,8,28,44,41 
1 .60-2.80 40,22,23,1 8,25,49, 10,2,27 
3.20-4.80 1,3,50,36,5,47,29,6 



0.20 
0.38 
0.68 



1.06 

2.13 
3.76 



0.49 
0.65 
0.60 



+0.03 
0.02 

+ 0.08 
0.05 

+0.07 
0.05 



0.07 

0.07 
0.07 



+0.02 
0.02 

+0.02 
0.02 

+0.02 
0.02 



0.52 
0.79 
0.69 



0.02 

0.01 
+0.03 

0.04 
+0.02 

0.03 



1.48 ±0.14 
2.17 ±0.28 
1.81 ±0.27 



Redshift-stacking for clusters with narrow temperature range: 1.20-1.34 keV 



z:0. 14-0.26 
z:0.28-0.30 



41,35,44 
51,28,8 



0.19 
0.29 



1.28 
1.27 



0.49 
0.56 



+0.03 
0.03 

+ 0.11 
0.07 



0.10 
0.14 



+0.03 
0.03 

+0.07 
0.05 



490+ 0012 

u.<*»u_ 00g 

52+ - 02 
u - oz -0.03 



1.49 ±0.18 
1.98 ±0.45 



Temperature-stacking for clusters with narrow redshift range: 0.23-0.33 



kT: 1.00- 1.60 
kT:1.70-2.80 



13,51,44,28,8,40 
22,23,25,18,10,27 



0.30 
0.30 



1.28 
2.10 



0.54 
0.55 



+0.07 
0.04 

+0.02 
0.01 



0.11 
0.07 



+0.05 
0.04 

+ 0.01 
0.01 



0.53 4 



3+0.02 
-0.02 

620+ 0013 
u.ozu_ 006 



1.83 ±0.32 
1.39 ±0.10 



Table 2. Results of the /3-modcl fit for the redshift-stacked and temperature-stacked clusters for all clusters and for redshift and temperature subsets. The 
f3 and r corc values in the fifth and sixth columns are for the /3-model with both r cor c and /3 freely-fitted. The values in the seventh column are for the 
/3-modcl with fixed r corc . All errors are la errors. 



in agree ment with studies of low re dshift groups ( s ee fo r 
example, iHelsdon & Ponmanl l200fj| and iMulchaev et~ai1 120031) . 
iHelsdon & Ponmanl ( 120001 ) attributed the fact that the slopes of 
group surface brightness profiles are flatter than the canonical slope 
(/3=0.67) for clusters, as a result of the effects of feedback from 
galactic winds on the intergalactic medium. The stacked data show 
a central excess, with f c — 1.56 ± 0.11, indicating that these sys- 
tems typically possess CCs. The individual profiles of z ~ 0.3 
clusters in Fig. [2] also support this results; the central excess factor, 
f c is greater than unity for 7 of the 12 clusters (25,51,22,8,28,40 
and 10) and an additional 4 clusters (44,13,18 and 23) have best fit 
f c > 1, but with error bars crossing unity. So, we conclude that 
58-92% of our systems at z ~ 0.3 have CCs. The fitted r CO rc for 
the stacked profiles is 0.06 x i? 5 oo- 



4.2 Trends of f c and f3 with redshift and temperature 

The main results of our analysis are the presence of trends in the 
value of /?, and in the incidence of cuspy cores, with redshift and 
temperature. Our CI clusters, as for any deep cluster survey, suffer 
from Malmquist selection effects, which result in increasing mean 
cluster luminosity with redshift, due to the fact that higher redshift 
cl usters are mo r e diffi cult to detect than nearby ones - see Fig. 3 
in lPacaud et al.l ( 120070 . Given the well-known correlation between 
X-ray luminosity and temperature, there is a corresponding ten- 
dency for more distant clusters in our sample to be hotter. This 
correletaion between z and T within our sample, makes it difficult 
to establish whether our observed trends in /3 and f c are evolution- 
ary effects, or whether they represent changes in cluster properties 
with system mass (and hence temperature). 

We attempt to address this issue in two ways: firstly, we can 
examine whether the trends we see (in both individual and stacked 
clusters) are stronger with respect to T or z. Secondly, we use the 
group of clusters at z ~ 0.3; subdividing these by temperature 
allows us to check for trends with T at essentially a single redshift. 



Similarly, we also extract a subset of our clusters which cover a 
rather narrow temperature range, but a larger spread in z. 

4.2.1 Trends of f c 

To investigate whether the trends we see in the redshift-stacked 
sets are affected by the T-z correlation in our sample, we select 
six clusters with similar temperatures (fcT=1.20 to 1.34 keV) but a 
relatively wide spread in redshift (z=0. 14 to 0.30). These are then 
divided into two subsets, each consisting of three systems: the first 
has 0.14 < z < 0.26 (average z=0.19) and mean temperature 
fcT=1.28 keV, the second has 0.28 < z < 0.30 (average z = 0.29) 
and kT= 1 .27 keV. The fit results for these subsets are shown in Fig. 
Hand in Tabled 

Similarly, we divided the twelve z ~ 0.3 clusters, which span 
a temperature range of 1.0-2.8 keV, into two temperature bins: 1.0- 
1.6 keV and 1.7-2.8 keV with six clusters in each. See Fig. [7] and 
the last two rows in Table|2]for the results of stacking these subsets. 
The 1.0-1.6 keV clusters have an average redshift z — 0.30 and 
average temperature fcT=1.28 keV, while the 1.7-2.8 keV clusters 
have the same average redshift z — 0.30 and kT=2.\ keV. 

The results from these subsets, shown in Fig. [6] and [7] rein- 
force the impression from Fig. [4] and \5\ that the increase in profile 
cuspiness is a function of z rather than T. In fact, the temperature- 
stacked subset (Fig.|7J actually shows a decline in f c with tempera- 
ture, whilst in the redshift-stacked subset it increases from 1.49 for 
the 0.14 < z < 0.26 clusters to 1.98 for the 0.28 < z < 0.30 
clusters (Fig.|6j. 

We also tested the f c -z behaviour in the individual profiles. 
In Fig. [8] we plot f c against redshift for the individual CI clusters. 
The figure shows that the high-z clusters tend to have larger-than- 
unity values of f c more often than clusters at lower redshift. We 
tested for a correlation in this plot, using the Pearson correlation 
coefficient, which has a value 0.40 for 26 points, corresponding 
to a Student t value of 2.12, which shows a poitive correlation at 
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Figure 4. X-ray surface brightness profiles of the redshift-stacked CI clus- 
ters with the associated la, 2a and 3a contours of the free r coro fit. The 
dashed lines are the fitted /3-model profiles with both r corc and freely 
fitted, while the solid lines are for the fitted profiles with free and r core 
fixed to 0.f05 x Rsoo- 



Figure 5. X-ray surface brightness profiles of the temperature-stacked C 1 
clusters with the associated la, 2a and 3a contours of the free r corc fit. The 
dashed lines are the fitted /3-model profiles with both r corc and freely 
fitted, while the solid lines are for the fitted profiles with free and r COTe 
fixed to 0.f05 x Rsoo- 
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over 95% significance (2-tailed test). To visualise the trend more 
clearly, we grouped adjacent data points into three bins and com- 
puted their weighted mean and the standard error. These binned 
results are shown as diamonds, though it should be noted that re- 
sulting values are sensitive to the choice of bin boundaries. 

In contrast, the / c -temperature plot in Fig. [9] shows no 
monotonic trend in f c with T, with a Pearson correlation co- 
efficient of 0.01. The binned values (diamonds) agree well the 
temperature- stacked results discussed earlier, where we noticed 
that intermediate-temperature clusters had the highest central ex- 
cess. 

The present work is the first study of the evolution 
of CCs within galaxy groups, although previous work (e.g., 
iHelsdon & PonmanluOOOl) has shown that CCs are common in X- 
ray bright groups at low redshift. Richer clusters have received 
much more study. CC clusters are found to be c ommon at low 
and moderate redshifts, see Bau er et al.l d2005f) . but lVikhlinin et al.l 
( 2007) found only a very small fraction of clusters at z >0.5 to have 
cuspy X-ray brightness profiles, w hich were taken as an indication 
of cool cores. ISantos et al.l (12008) found moderate CC clusters out 
to 2 = 1.4, but noted an absence of strong CCs at redshifts higher 
than 0.7. 

Our results therefore suggest that groups behave differently to 
clusters, in that cuspy cores are actually more prominent at higher 
z in these poorer systems. How can we understand this differ- 
ence? One possibility is that the central excess seen in groups at 
moderate-high z is not due to CCs at all, but to the presence of 
central AGN. We can immediately rule out the possibility that the 
effect is due to just a few groups with bright AGN contaminating 
our stacked profiles by noting (cf. Fig. [8) that a central excess is 
seen in the majority of systems at z >0.3. Hence, any effect from 
central AGN would have to be widespread and moderate. 

There would be significant spectral differences between cen- 
tral excesses generated by CCs and AGN, since the thermal emis- 
sion from cool cores is much softer than the X-ray spectra of AGN. 
The limited statistics for individual clusters in our sample do not 
permit us to investigate whether the core emission is soft or hard. 
However, this can be investigated using the stacked data. We there- 
fore repeated the stacking analysis for intermediate- and high- 2 
clusters using X-ray images derived from the hard energy band, 
2.0-4.5 keV. The results are shown in Fig.[l0] Comparing with the 
corresponding soft (0.5-2.0 keV) band profiles in Fig. [4] we no- 
tice the disappearance of the central excess above the /3-model in 
the hard-energy profiles in both the intermediate and high redshift 
stacks. This provides strong evidence that this central excess does 
not arise from AGN in cluster cores. 

Assuming that the cuspy profiles really do indicate the pres- 
ence of CCs, the observation of such cores in groups at high redshift 
can help to constrain the reasons for their absence in high z clusters. 
The decline in CC clusters with redshift could result from disrup- 
tion of CCs due to th e h igher merging rates at high redshifts (e.g., 
ICohn & Whitej|20 05 and Jeltema et alj|2005l) or from the effects of 
preheating ( McCarthy et al. 2004), which can raise the gas entropy 
and prevent cooling. In the latter case, the impact of a given entrop y 
boost will be larger in cooler systems (see e.g. lBorgani et al. 2005), 
especially at high redshift, so one would expect if anything to see 
a drop in the incidence of CCs within groups at high 2, at least as 
large as pronounced as that in richer clusters. 

The hypothesis that CCs are destroyed by cluster mergers ap- 
pears more promising, since t his effect might be stronger in more 
massive systems. For example. lBurns et alj (12008) find that CCs are 
more common in low mass clusters, and attribute the lack of CCs in 



more massive systems to their destruction by early mergers in sys- 
tems destined to grow into large clusters. On the other hand, these 
authors caution that their model does not reproduce the observed 
reduction with redshift in CCs within massive clusters. In fact, no 
numerical simulations have yet succeeded in adequately reproduc- 
ing the properties of cluster cores. 



4.2.2 Trends of (3 

Whilst it seems quite clear, as discussed above, that the trend in 
central cuspiness is primarily related to redshift, rather than temper- 
ature, the situation with regard to /3 is not so straightforward. In the 
stacked datasets, comparison of Fig. [4] and \5\ suggest that the rela- 
tionship with redshift is stronger: /3 rises monotonically through the 
three redshift intervals, whilst in temperature the only clear result 
is that the cool systems have lower /3. In contrast, the fits to individ- 
ual clusters ("Fig. II llandll2b show a more pronounced trend with 
temperature than with redshift. The stacked subsets of the narrow 
temperature and redshift ranges (Fig.[7]and[6} produce ambiguous 
results: /3 increases with temperature if r corc is fixed, but with z if 
r corc is left free to fit. 

So, on the basis of our data alone, we are unable to say whether 
the general trend in /3 is driven by the temperature or the redshift. 
However, evidence from studies of low redshift groups and clus- 
ters is very relevant here. Such studies provide clear evidence of 
a positive correl ation between j3 and temp erature in local systems 



ind tempt 
n (2004) 



systems - e.g., Osmo nd & Ponmanl d2004l) and Fig. 7 in the study 
of Croston et al. (2008), who analysed clusters with redshift < 0.2. 
Combining these previous results with ours, favours a trend in /3 
with temperature (and hence cluster mass) rather than an evolu- 
tionary effect. Our results then demonstrate that this trend is still 
present in groups and clusters at z ~ 0.3. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we used XMM-Newton observations of 27 X-ray se- 
lected galaxy clusters spanning the redshift range (0.05 < z < 
1.05) to study the spatial properties of their ICM. Most of these 
clusters fall in the realm of low-mass clusters or groups, with 
ICM temperatures from 0.63 to 4.80 keV. The XMM data provide 
typically a few hundreds X-ray source counts. We extracted and 
vignetting-corrected the profiles to 3xR,5oo where they flattened 
and reached the photon background values, which were estimated 
locally for each cluster. 

In addition to the individual profiles, we also stacked the pro- 
files into three redshift and temperature bins. To explore the effects 
of Malmquist bias, we further stacked clusters with similar red- 
shifts/temperatures into two subsets each with different averaged 
temperature/redshift. Both individual and stacked profiles were fit- 
ted with blurred (to account for the PSF errors of the XMM-Newton 
cameras) /3-models with both free and fixed core radii. The fixed- 
?~corc /3-model fits were used to test whether a profile showed evi- 
dence of a cuspy core, making this study the first to probe the evo- 
lution of CCs out to 2 > 0.3 within poor clusters. 

Our main conclusions are: 

• We find that 54% of our sample show evidence for cool cores, 
in the form of a central excess (at > la significance) above a stan- 
dard /3-model. 

• For the free-r corc fits to individual clusters, the median value 
of f3 is 0.61, and the median r coro is 0.08xRsoo- 
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Figure 6. X-ray surface brightness profiles of stacked C 1 clusters with nar- 
row temperature range (1.20-1.34 keV), grouped into two redshift bins: 
0.14 < 2 < 0.26 (top panel) and 0.28 < z < 0.30 (middle panel). 
The bottom panel is the lcr, 2cr and 3a contours. The dashed lines are the 
fitted /3-model profiles with both r core and /3 freely fitted, while the solid 
lines are for the fitted profiles with free (5 and r CO rc fixed to 0. 105 X R500 ■ 



Figure 7. X-ray surface brightness profiles of stacked C 1 clusters with nar- 
row redshift range (0.23-0.33), grouped into two temperature bins: 1.00- 
1.60 keV (top panel) and 1.70-2.80 keV (middle panel). The bottom panel 
is the lcr, 2<r and 3<r contours. The dashed lines are the fitted /3-model pro- 
files with both r corc and f3 freely fitted, while the solid lines are for the 
fitted profiles with free /3 and r CO rc fixed to 0.105 X Rsoo- 
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Figure 8. Central excess factor, f c (of the fixed r core fits) plotted against 
redshift of the individual CI clusters. f c is defined as the ratio of the ob- 
served surface brightness to the predicted (model) surface brightness within 
O.O2XR500 (first radial bin). A value of f c above 1 is an indication of 
a cool core cluster and vice versa. The positions and sizes of the diamonds 
represent the weighted means and the standard errors of the weighted means 
of the points as described in text and appendix lAl 
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Figure 9. Central excess factor, f c (of the fixed r corc fits) plotted against 
temperature of the individual CI clusters. The positions and sizes of the 
diamonds are calculated as described in appendix lAl 



• For the fixed-r corc fits to individual clusters, the median ft is 
0.63. 

• Twelve systems in our sample (with kT — 1.69 keV) have 
z ~ 0.3, allowing us to characterise the X-ray surface brightness 
profiles of intermediate redshift X-ray selected groups. The free-fit 
parameters to the stacked data from these 12 systems gives /3 = 
0.51 ± 0.01 and r core =0.06 ± 0.01xR 50 o- This stacked profile 
indicates the presence of CCs (/ c = 1.56 ± 0.11), with 7 of the 
12 systems showing a significant central excess in their individual 
profiles. 

• Stacked and individual profiles for our sample of poor galaxy 
clusters show that the CCs do not disappear at high redshift, but 
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Figure 10. Stacked profiles for the intermediate- and high-redshift cluster 
subsamples extracted from hard band (2.0-4.5 keV). These can be compared 
directly with the corresponding soft band stacks shown in Fig. [4] and shown 
that the central excess above the fixed-r corG /3-model is not present in the 
hard band. 



rather become more prominent, though one would like to confirm 
this result with higher spatial resolution observations. 

• The slope parameter, /?, shows a positive trend with both red- 
shift and temperature in our data. Combining this results with pre- 
vious findings, we incline towards a trend with temperature (and 
hence mass) rather than redshift. The present study then demon- 
strates for the first time, that the f3-T trend seen at low z is also 
present in groups and clusters at z > 0.3. 
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Figure 11. /3 values (of the free r corc fits) versus redshift of the individ- 
ual CI clusters. The positions and sizes of the diamonds are calculated as 
described in appendix lAl 
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APPENDIX A: STANDARD ERROR ON WEIGHTED 
MEAN IN PRESENCE OF REAL SCATTER 

When averaging data of variable statistical quality, a more robust 
mean is obtained by weighting the averaged values by their in- 
verse variances. Standard formulae for the standard error on such 
a weighted mean assume that statistical errors represent the only 
source of variance. For our application, this is not true, since 
there are real cluster-to-cluster variations, in additional to statist- 
cal scatter. Here we derive an expression for the standard error of a 
weighted mean in these circumstances. 

For a data set x% — Xi, X2, •-, x n with variable statistical er- 
rors o~i — err, (72, .., a n , the weighted mean x is 



and finally the standard error of the weighted mean (SEWM), SE S , 
is 



SE a = 
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This paper has been typeset from a TgX/ KTgX file prepared by the 
author. 
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where Wi = 1/a 2 are the weights. This weighted mean will prop- 
erly take into account the varying statistical quality of the data. 
The variance in this the weighted mean, var(x), is 



var(x) = 



N w i var(xi) 

i=l 



(A2) 



E* 



In the presence of real, non-statistical scatter in the x values, 
the expected variance for the ith data point is 

var(xi) = ((xi - fi) 2 ) - af + at, (A3) 

where a t is the true (non-statistical) variance of the population and 



((Xi - m) 2 ) = t ((Xi - xf) . 

So an estimate a 2 of a 2 is obtained from 
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and the variance of the weighted mean becomes 
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Substituting for Wi = I /a 2 we get 
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